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We study the attractive interactions between rod-like charged polymers in solution that appear 
in the presence of multi-valence counterions. The counterions condensed to the rods exhibit both 
a strong transversal polarization and a longitudinal crystalline arrangement. At short distances 
between the rods, the fraction of condensed counterions increases, and the majority of these occupy 
the region between the rods, where they minimize their repulsive interactions by arranging them- 
selves into packing structures. The attractive interaction is strongest for multivalent counterions. 
Our model takes into account the hard-core volume of the condensed counterions and their angu- 
lar distribution around the rods. The hard core constraint strongly suppresses longitudinal charge 
fluctuations. 

PACS numbers: 61.20.Qg, 61.25.Hq, 87.15.Aa 
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Strongly charged polymers precipitate from a dilute so- 
lution into compact structures when high-valence coun- 
terions (oppositely charged particles) are added to the 
solution jlj-g]. The counterions experience strong elec- 
trostatic attractions to the backbone of the chains, and a 
finite fraction of them "condense", i.e., are found within 
short distance from the chains |Q . Counterions are more 
attracted to compact chains or aggregates of rod-like 
chains. This creates the possibility of a transition from 
single chains with a small number of condensed coun- 
terions to almost neutral aggregates of chains, or even 
mono-molecular collapse in the case of flexible polymers. 
These aggregates are stable only when the internal ar- 
rangement of the counterions within them, provides a 
strong enough cohesive energy. 

In this letter we study the attraction between two rod- 
like polyelectrolytes. We show that it is essential to in- 
clude the size and angular degrees of freedom (around the 
rods) of the counterions as well as the discrete nature of 
the charge along the polyelectrolytes to find the origin 
and strength of the counterion mediated attraction. Our 
work suggests that these factors are also crucial in de- 
termining the collapse of flexible and semi-flexible poly- 
electrolytes recently studied in Refs.[8-12]. Experimental 
observations show that the size of the precipitating par- 
ticles is indeed a relevant parameter in the problem . 

It has been argued that longitudinal charge fluc- 
tuations resulting from the thermal motion of 
point counterions induce attractions between rod-like 
polyelelctrolytes jl3],[L4| and induces "buckling" of semi- 
flexible polyelelctrolytes |ll|,|l2|. Here, we show that 
such charge fluctuation are suppressed when the hard 
core volume of monomers and counterions are taken into 
account. Instead, we find that the counterions arrange- 
ment around the rods create a non-zero transversal po- 
larization as the distance between the chains decreases. 
At very short distances between the rods we find strong 
longitudinal correlations but only at very short wave- 
lengths, implying a crystalline state along the rod, re- 



inforcing, for the case of multivalent counterions, the 
attractive interactions due to polarization. 

The crystalline structure of the counterions when the 
rods are at short distance from each other has been ob- 
served in simulations by Gr0nbech- Jensen et al. [jig] , and 
has been theoretically proposed by Arenzon et al. 
and Shklovskii jl7| . These previous theoretical works re- 
tain some of the small size effects of a realistic system, 
but again use the assumption of negligible size counteri- 
ons. These models do not reveal the polarization effects 
that appear when the angular degrees of freedom around 
the rod are considered. 
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FIG. 1. The diameter of the rod, the size of the coun- 
terions and the basic spacing between charges all have 
magnitude b. We use two representative states for each 
monomer in each rod. All occupation is assumed to take 
place parallel to the centers of the monomers, so that the 
possible places for occupation form the lattices 1, 2, 3, and 
4, here shown from the top. 

The polymer chains are modeled as rigid rods formed 
by TV repeat units, each with charge +1 centered at each 
of the monomers, as shown in Fig. 1. We label the pair 
of rods as A and B. The diameter of the rods and the 
spacing between charges are both b. The total length of 
the chains is L — Nb. We will only consider the case in 
which the chains are parallel and separated by a distance 
r measured from center to center. The counterions carry 
a charge —Z, and for simplicity we take them to be of 
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diameter b as well. Overall charge neutrality of the sys- 
tem implies that there are N/Z counterfoils per chain. 
We assume a concentration of chains c, to be such that 
the average distance between rods (if non-interacting) is 
also L, so that c = 1/L 3 . 

A layer of condensed counterions surrounds each chain. 
The location of a particular counterion is given by its po- 
sition along the rod z, a radial distance r from the rod, 
and its angular position 0. We restrict the z coordinate 
of the condensed counterions to take values coinciding 
with those of the centers of the monomers. The radial 
distance form the rod for all condensed counterions is 
assumed to take always the value r = b. The angular 
variable is very important as it carries information about 
the local polarization of the rod-counterions system. A 
suitable simplification that retains this information con- 
sists on collapsing the range of angular positions to only 
two specific locations. These positions are located at an- 
gles of 7r/4 and —3tt/4 with respect to the plane that 
contains both rods, and are labeled as shown in Fig. 1. 
These ad-hoc positions are meant to represent locations 
of counterions on one rod that look towards or away from 
the other rod, and have the convenience of allowing us 
to consider distances between the centers of the rods as 
small as y/2b without worrying about overlapping of the 
counterions that would require explicit introduction of 
hard-core repulsions. The choice of the angles does not 
significantly change the results. In short, we have mod- 
eled the condensed layers as four linear lattices parallel to 
the rods, whose sites can be occupied by the counterions. 

Since our selected geometry takes care of the hard-core 
interactions between the particles of the system, we can 
construct a Hamiltonian with only electrostatic interac- 
tions. This is given by: 

We measure energies in units of fcgT, and the prefactor 
Ib is the dimensionless ratio of the Bjerrum length to 
the monomer size e 2 / ebksT , with e being the electron 
charge, k B the Boltzmann constant T the temperature 
and e the dielectric constant of water. <f> is the electro- 
static potential created by the charged rods. Ignoring end 
effects from the rods, this potential is well approximated 
by 4>{t) = 2l B (\n(L/r A ) + \n(L/r B )), where r A ,r B , are 
the distances from the point r to the axis of the rods A 
and B, respectively. For the condensed counterions, it is 
better to change to a local charge representation. Each 
site of the four lattices can be occupied by one counterion 
(we neglect multiple occupation), and thus it will carry 
a charge qi{n) that can be either — Z or zero. The index 
i is the is the (longitudinal) position in the lattice, that 
can range from 1 to N. 

Since the number of monomers is large, we expect that 
the total number of condensed counterions for a given 



inter-rod separation will have a narrowly peaked proba- 
bility distribution. We construct a free energy that as- 
sumes that the number of counterions condensed to each 
of the lattices is fixed, and then we will find the minimum 
with respect to the occupation numbers. Lattice i carries 
a fraction /j of the number of counterions per rod N/Z, 
and the fraction condensed to rod A is Ja = fi + f2, etc. 

The free counterions form a dilute charged gas that oc- 
cupies a volume V = 2L(L 2 — 4b 2 ). Because of their low 
density their contribution to the free energy from corre- 
lations and screening is negligible. A free counterion at 
a distance r c from the center of the system feels a poten- 
tial given approximately by <j> = 2(lg/e)(2 — Ja~ /b)(1 — 
(r 2 /2L 2 )) ln(L/r c ), that arises from the effective (uncom- 
pensated) charge of the two rods and a cylindrical shell 
of uniformly distributed free counter-ions. Averaging this 
potential over the volume, and adding the entropic con- 
tribution, we obtain the free energy per monomer due to 
the free counterions: 

F f = -^(2 -j A - Ib) (M(2 -f A ~ f B )Nb 3 /V) - l) 

+ \^-fA-f B ) 2 ^-\lB^ (2) 

To obtain the contribution from the condensed coun- 
terions, we consider first a high-temperature approach 
in which we add fluctuations to a uniformly distributed 
state. Given a condensed fraction /j at the lattice i, 
there are fiN/Z occupied sites. The statistical sum over 
all states satisfying this restriction can be replaced at 
high temperatures by Gaussian integrations over a set 
of continuous local charge variables. The charge qi(n) is 
represented by a density p%{n) with mean — /j and vari- 
ance of = Z 2 (fi/Z)(l — fi/Z). This form of the variance 
is consistent with our assumption of a maximum of one 
counterion per lattice site, as required by our geometry 
and hard-core constraints . We can pass to a discrete 
Fourier representation of the local charge of the form 

Pi(n) = - fi + E pi{k) exp(ikn) (3) 

where the only Fourier modes considered are of the form 
k = ±2wm/N, with m ranging from —N/2 to N/2. This 
transformation diagonalizes the part of the Hamiltonian, 
Eq. (I), that corresponds to the condensed counterions, 
and leads to the free energy per monomer: 

f c = \ E fiKfi + i^F E lndct P + sv ( fc )i 

i,j fc/0 

The first term is the contribution from the zero modes 
where the interaction matrix between lattices is Vf\ = 
21b ln(L/rij) with the distance between the axis of 
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the lattices, and the diagonal terms are given by = 
2Zb m(£/&). The second term is the result of the Gaus- 
sian integration over fluctuations, where I is the identity 
matrix, S is a diagonal matrix with entries cr|, and the 
interaction terms for the Fourier modes are of the form 
Vij(k) = 2lBK (krij), and Vu(k) = —2lBci(kb), with K 
the modified Bessel function, and ci the cosine integral 
function. The integrals that define the diagonal expres- 
sions are evaluated using the lattice spacing b as short 
distance cutoff so that there is no need for self-energy 
subtractions. The finite sum over modes also corresponds 
to consider distances larger than one lattice spacing. The 
last term is the entropic contribution of placing fcN/Z 
counterions in a volume V c = Lb 2 . 

A posteriori we have found that at short distances 
between the rods the condensed counterions are almost 
fully polarized, occupying the states that face towards the 
other rod. The fraction of condensed charge approaches 
1 in these sites, while the outward states are almost com- 
pletely depleted. We can look then in more detail at the 
calculation of the free energy for the case of almost full 
occupation of one of the lattices / — > 1, so that the vari- 
ance of the fluctuations becomes a 2 » (Z — 1). The 
diagonal interaction term for the highest Fourier modes 
k w ±7r/& is negative and of order Is, corresponding to a 
decrease in energy from vacating a state next to an occu- 
pied one. Since typical values for Ib at room temperature 
are always larger than 1, these diagonal terms dominate 
the interaction matrix V, and therefore the matrix I+SV 
has diagonal elements 1 — (Z — 1 )Ib- When the condition 

(Z-1)1 B >1 (5) 

is satisfied, the matrix acquires negative eigenval- 
ues, making the determinant divergent and the high- 
temperature approximation incorrect. Clearly, this oc- 
curs in most cases, except for mono-valent counterions 
(Z = 1) or for very weakly charged polyelectrolytes for 
which Ib <C 1. The divergence in the determinant for 
multi-valent counterions signals the onset of crystalliza- 
tion, and thus the free energy should be calculated on 
the basis of a dominant crystalline ground state. Once 
we use the correct ground state for both mono- and multi- 
valence cases, it can be shown that the corrections from 
fluctuations are very small, that is, the contribution of 
the determinant in Eq. 4 becomes negligible once the di- 
vergent modes arc subtracted. 

The proper ground state for each of the inner lattices, 
when approaching full condensation, is clearly given by 
an arrangement in which counterions are placed one in 
every Z sites. (This is the case already for Z = 1.) As 
the rods approach each other, the inner lattices interact 
strongly, but they do not destroy the ground state ar- 
rangement. Instead, they can choose a location of the 
occupied sites, so as to minimize their repulsive interac- 
tion. For Z = 2, for example, one expects one of the 



internal lattices to be filled in the even sites, while the 
second in the odd ones. 

The calculation of the free energy for the multi- valence 
case can be carried out using the results for the mono- 
valent case which are given by Eq. (4). We simply renor- 
malize the lattice spacing to Zb and reduce the avail- 
able sites by a factor of Z. The elements of the inter- 
action matrix for the ground states are now given by 
V-j — 21b ln(L/r^ -), with a modified distance between 
the lattices = (r?- + {bZ/2) 2 ) 1 / 2 that takes into ac- 
count the mismatch between the occupied sites in the lat- 
tices. The diagonal term is simply V™ = 21b ln(L/Zb). 

In Fig. 2 we present the calculated free energy of the 
system as a function of the distance between rods. The 
numerical values for the constants of the system are 
T = 300, b = 1.8A°, so that l B = 4.1, and N = 10 5 . 
There is a well of attraction for the mono-valent case of 
about .5kT, which is not enough to bind the rods, and 
further, the local minimum at short distances turns out 
to be of higher energy than the self-energy of two rods 
separated by a distance L. For Z = 2 and 3, the depth 
of the well is of the order of ksT, and the energy there 
is lower than their respective reference states. 
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FIG. 2. Results of the minimization of the free energy, 
for parameters such that Ib = 4.f, for valences Z = 1,2, 
and 3. Energies are shown in units of fcsT and distances 
in units of b. The zero of the energy in each case is chosen 
to match the energy at large separations. 

Fig. 3 presents the total amount of condensed counteri- 
ons for one of the rods, j a = fi + fi- This is always very 
close to 1 for short distances and reaches a value near 
the Manning limit [Q, /= (1 — 1/ZIb) at large separa- 
tions. We measure the overall polarization p of the rods 
by the ratio of the difference between the occupation of 
the inward and outward positions to the total amount of 
condensed charges, thus for rod A: 



Results for the amount of polarization are shown in 
Fig. 4. It is clear that the charge will be perfectly bal- 
anced when the presence of the second rod is not felt. 
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On the other hand, when the chains are in close contact 
is natural for the counterions to occupy the inner sites to 
be able to interact with the positive charges of both rods, 
even at the expense of interacting with other condensed 
counterions. What it is surprising is that both the polar- 
ization and the extra condensation do not decay quickly, 
and it is necessary to set the distance between the rods 
to its maximum value L to recover the Manning limit 
and a symmetric state. A good test of the validity of 
this theory will be the measurement in simulations and 
experiments of the transversal polarization of the rods. 



is an important attraction between the rods when the 
counterions are multi-valent. We found that the finite 
size of the counterions and their angular degrees of free- 
dom are essential to determine the nature and strenght of 
the counterion mediated attractions in rigid-rods, and we 
expect this to be also the case in flexible and semiflexible 
polyelectrolytes. 
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FIG. 3. The total fraction of condensed counterions to 
each of the rods, as a function of the separation between 
rods for valences Z = 1, 2, and 3. The scale in which the 
condensed fraction decays to a single rod value is of the 
order of the size of the rod L = 10 b. Note the change of 
scale with respect to Fig. 2. 
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FIG. 4. Polarization of the rods as a function of the 
separation between rods for valences Z = 1,2, and 3. At 
very short distances, the polarization is almost complete. 
Away from the region of strong attraction, the polariza- 
tion still remains important and will decay to zero only 
at distances of order L. The plots for Z = 2,3, overlap 
almost completely. 

In summary, we have shown that the interaction be- 
tween two charged rod-like polymers generates a strong 
transversal polarization of their condensed charges and 
that at short distances the two rods are strongly driven 
towards higher counterion condensation. This forces the 
counterions to crystallize and then to organize their re- 
spective crystals into a packing structure. The final result 
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Our variance differs from the expression used in models 
that assume point-size counterions, refs. 0,0, where 
the variance is proportional to /j. A variance proportional 
to fi implies, at large condensation fractions, frequent 
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multiple occupation at a monomer site that is inconsis- between the rods, 

tent with the reduced space available to the counterions 
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